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Abstract. The mechanism of corc-coUapse supernova explosions must draw on 
the energy provided by gravitational collapse and transfer the necessary fraction 
to the kinetic and internal energy of the ejecta. Despite many decades of concerted 
theoretical effort, the detailed mechanism of core-collapse supernova explosions 
is still unknown, but indications are strong that multi-D processes lie at its 
heart. This opens up the possibility of probing the supernova mechanism with 
gravitational waves, carrying direct dynamical information from the supernova 
engine deep inside a dying massive star. I present a concise overview of the 
physics and primary multi-D dynamics in neutrino-driven, magnetorotational, and 
acoustically-driven core-collapse supernova explosion scenarios. Discussing and 
contrasting estimates for the gravitational-wave emission characteristics of these 
mechanisms, I argue that their gravitational- wave signatures are clearly distinct 
and that the observation (or non-observation) of gravitational waves from a nearby 
core-collapse event could put strong constraints on the supernova mechanism. 

PACS numbers: 97.60.Bw, 97.60.Jd, 97.60.-s, 97.10. Kc, 04.30. Db, 04.40.Dg 
Submitted to the GWDAW13 special issue of CQG on May 17, 2009. 

1. Introduction 

In their 1934 paper Baade and Zwicky suggested that supernovae (SNe) are 
caused by the release of gravitational energy in the "collapse of ordinary stars to 
neutron stars." Now, seventy-five years after this initial proposition and after much 
observational and theoretical effort, Baadc's and Zwicky's statement still holds and is 
the fundamental paradigm of corc-coUapsc SN theory. 

Massive stars (~8 — 10 Mq ^ M < 100 M0) live fast nuclear-burning lives, 
consuming their nuclear fuel in a cosmic blink of only '^100 million — 1 million years 
(r/r© « {M/Mq)~'^-^). When nuclear burning ceases in the core, a massive star has 
evolved to a red or blue supergiant (radius ~ few x 10^'^ cm or ~ 0.1 — 1 x 10^^ cm, 
respectively) and is comprised of a central, compact Chandrasekhar-mass electron- 
degenerate iron core (R ~ 1500 km) that is surrounded by an onion-skin structure 
of shells of material with successively lower mean atomic weight. Supported only 



Probing the SN Mechanism with Gravitational Waves 



2 



by increasingly relativistically-dcgcncratc electrons, the iron core eventually exceeds 
its effective Chandrasekhar mass. Catastrophic gravitational collapse sets in and is 
accelerated by electron capture on nuclei and free protons as well as photodissociation 
of iron-group nuclei. In a matter of a few hundred milliseconds, the inner core is 
compressed from densities of 10^ — 10^" g cm~^ to nuclear density (~ 2.6 x 10^"' g cm~'^) 
where the nuclear equation of state (EOS) stiffens and stabilizes the inner core. The 
latter overshoots its new equilibrium, then rebounds into the still collapsing outer 
core. This core bounce results in the formation of a hydrodynamic shock at the 
interface of inner and outer core. The shock propagates outward in radius and mass 
through the supersonically infalling outer core material, leaving behind the unshocked 
~ 0.5 — 0.7 Mq inner core which forms the core of the hot protoneutron star (PNS). 
The PNS rapidly gains mass through accretion (M is of the order of 1MqS~^) and 
radiates away its gravitational energy (~ 300 B for a 1.4 Mq NS; 1 [Bjethe = 10^^ erg) 
in neutrinos as it evolves to a cold NS on a timescale of tens of seconds. 

In early theory and SN models, a prompt hydrodynamic explosion would follow shock 
creation Q , removing the stellar envelope and providing the fantastic optical display 
of a SN with total kinetic, internal and electromagnetic energy of ~ 1 Ep. 

More advanced theory and simulations with more complete and accurate treatment 
of the EOS, neutrino physics and transport show that the prompt shock 
fails: Photodissociation of infalling heavy nuclei into nucleons (at a cost of ~ 
8.8 MeVbaryon^^ for iron-group nuclei) and neutrinos that stream away from the 
optically-thin postshock region sap the shock's energy. The shock slows down quickly 
after its launch, stalls, and turns into an accretion shock at a radius of 100 — 200 km. 
For a SN explosion to occur, there must operate a mechanism that, in some fashion, 
transfers and deposits a fraction of the tremendous gravitational energy of collapse into 
the immediate postshock region, reviving the shock and endowing it with sufficient 
kinetic energy to make a core-collapse SN. 

What is the core-collapse SN mechanism? This is the fundamental question and 
primary unsolved problem of core-collapse SN theory. Finding observational evidence 
for the SN mechanism by classical astronomical means from radio to 7 wavelengths 
is difficult, since all pre-explosion dynamics takes place deep inside the stellar core, 
completely shrouded from view in the electromagnetic (EM) spectrum. Hence, EM 
observations yield only secondary observables, such as progenitor type and mass, 
explosion morphology and energy, ejecta composition, neutron star properties and 
birth kick. Primary, direct "live" information from deep inside the supernova engine 
is carried only by neutrinos and gravitational waves (GWs) which both can propagate 
from their emission sites to observers on Earth virtually without interaction with 
intervening matter. 

GWs have not yet been observed directly, but an international network of 
interferometric detectors (LIGO, VIRGO, GEO, TAMA) is taking data with sufficient 
sensitivity to observe GWs from a nearby galactic core- collapse SN (e.g., 0, H). 
Neutrinos, on the other hand, were detected from SN 1987A, confirming the basic 
theory of stellar collapse and neutrino emission. However, in part due to the small 
number of detected neutrinos, this observation left unanswered the question of the SN 
mechanism. 

In this contribution to the proceedings of the 13th Gravitational- Wave Data Analysis 

X Inferred core-coUapsc SN energies fall in the range of 0.1 — 10 B Q, but energies around IB are 
most frequently seen. 
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Workshop (GWDAWIB), I delineate the current set of proposed core-cohapse SN 
mechanisms. I discuss the underlying physics of each mechanism and limit my 
focus to the dominant multi-D processes active and leading to GW emission in a 
given mechanism. Furthermore, I argue that the GW signatures of the various SN 
mechanisms are distinct and that the observation or non-observation of GWs by 
current and/or future advanced or third-generation GW observatories from a nearby 
SN can put significant constraints on the explosion mechanism. This idea was put forth 
first in [J] and is elaborated here. In section [21 I discuss the neutrino mechanism and 
its GW signature. This mechanism is currently favored as the explanation for core- 
collapse SN explosions in garden- variety slowly or nonrotating massive stars. Section[3] 
is devoted to the magnetorotational mechanism and its GW emission. This mechanism 
may be active in rapidly rotating massive stars. In section [U I discuss the recently 
proposed acoustic mechanism of Burrows et al. 0, [3| which relies on the excitation 
of PNS pulsations and their damping via acoustic waves that steepen to shocks and 
efficiently deposit energy in the postshock region. In section [51 I summarize and draw 
conclusions. 



2. The Neutrino Mechanism 

The neutrino- driven mechanism (or simply, neutrino mechanism), is founded on the 
notion that if only a small fraction of the iJgrav ^ 300 B of gravitational energy 
liberated in collapse and emitted in neutrinos and antineutrinos of all flavors was 
deposited in the stellar mantle, it could be driven to explosion. The neutrino 
mechanism, based on the charged-current absorption of Vg and on neutrons and 
protons, respectively, was first proposed by Colgate & White Q and Arnett Q and 
was initially envisioned to be efficient and direct, that is, expected to drive the shock 
out before it has time to stall. However, simulations with improved numerics and 
input microphysics showed that neutrino energy deposition ( "neutrino heating" ) is 
actually rather inefficient and the shock always stalls first and may be revived only on 
a timescale of hundreds of milliseconds in the delayed neutrino mechanism proposed 
by Bethe & Wilson The amount of energy to be delivered by the mechanism, 

-E-nicch, is set by the sum of asymptotic explosion energy of ~ 1 B and the gravitational 
binding energy of the mantle (also of the order of 1 B) that must first be overcome. 
Hence, £'mcch/£'grav ~ 1% is the fraction of the total energy that must be transferred 
to the mantle by neutrinos to power an explosion. This 1% must be provided to the 
shock within a short, < 1 — 2s time frame after bounce in order for the onset of 
explosion to occur before accretion has pushed the PNS over its maximum mass and 
for the SN to leave behind a NS mass consistent with observations. Typical Ve and 
luminosities in the postbouncc phase are L^^ « Lp^ w 10Bs~^ and, hence, assuming 
that the explosion must start within ~ 1 s, the necessary heating efficiency is of the 
order of - 10%. 

As compelling and promising the neutrino mechanism may seem, the most advanced 
simulations predict its failure in ordinary massive stars with masses > 10 Mq when 
spherical symmetry (ID) is assumed and fluid motions are restricted to the radial 



direction [llj, [ij, [l3| . Only the lowest-mass massive stars with 7 Mq ^ M < 10 Mq (S- 



AGB stars that end their nuclear burning lives with 0-Ne cores) may explode, though 



sub-energetically, by the ID neutrino mechanism alone [IJ, [15|- In multi-D, neutrino- 
driven convection [13, in combination with the standing-accretion-shock 
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Figure 1. Postbounce supernova cores of the nonrotating model s20.nr (left 
panel) and the rapidly-spinning model s20.7r (right panel) of the angle-dependent 
neutrino-radiation-hydrodynamics simulations of Ott et al. l25f at ~ 215 ms after 
bounce. Shown are the specific entropy distributions and velocity vectors are 
superposed to visualize flow patterns. Key aspects of the postbounce situation 
are highlighted in model s20.nr in the left panel. 



instability (SASI) (e.g., [1^, 20, 21 1 and references therein) can enhance the heating 



efficiency [2^|23|. Yet, the presently most complete axisymmetric (2D) simulations 
still fall short of producing robust 1-B neutrino-driven explosions @, [22, [24, 2^ and 



hope for the neutrino mechanisms rests on the additional degrees of freedom provided 
by 3D dynamics to facilitate robust explosions in nature and in future detailed 3D 
simulations [2a, 13, 2^ 

Figure [T] displays color maps of the specific entropy distribution in snapshots taken 
at ~ 215 ms after bounce of the simulations of Ott et al. [2^ that started with a 
20-Mq presupernova model and were carried out without rotation (model s20.nr) and 
with rapid precollapse rotation (model s20.7r, Qq = Trrads"^). Simulation details are 
explained at length in (25| . The nonrotating model is depicted in the left panel and 
the various features of the postbounce situation are highlighted by text and arrows. 
Neutrinos are trapped in the dense PNS core, but leak out from the neutrinospher^i 
at the PNS surface. Neutrino cooling {Q~ oc T^, e.g., [131) dominates over neutrino 
heating (Q^ oc L^r~'^{e^), interior to the gain radius of 80 — 100km denoting 
the location beyond which (5+ > Q^. The net neutrino heating is strongest near the 
inner boundary of the gain region and decreases outwards. Hence, a negative entropy 
gradient develops, resulting in instability to overturn and convective eddies in the 
gain region are clearly visible. Also discernible is the deformation of the shock front 
due to the early phase of the SASI which drives the growth of low- mode {£ = 1,2) 
perturbations in an advective-acoustic cycle of laterally propagating sound waves and 
advection of pressure perturbations (e.g., (2ll | and references therein). 

§ The neutrinosphere is defined as the location at which the neutrino optical depth is 2/3. The 
radius R^, of the neutrino sphere is strongly dependent on neutrino energy. 
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2.1. Convection and SASI in Nonrotating and Rotating SN Cores 



In the core-collapse SN context, convective instability can generally be driven by 
negative gradients in entropy s or lepton fraction Yi ~ Ye + Y^. Rotation or more 
specifically, a positive specific angular momentum gradient, can stabilize the flow 
against overturn. This is expressed in the Solbcrg-H0iland condition for instability. 
In the equatorial plane and neglecting lateral gradients, this condition is given by 

1 d 



dr 



■f < 



(1) 



Here, j = flr^ is the specific ang ular momentum and iV^ is the Brunt- Vaisala 
frequency, in the SN context [H, l29t 
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where g is the gravitational acceleration, 7 = dlnp/dln p\s, and dp/ds and dp/dYi 
are generally positive. In the right panel of figure [1] the entropy distribution and the 
superposed velocity vector field of the rapidly spinning (PNS period ~ 2 ms) oblate 
postbounce SN core of model s20.7r is shown at ~ 215 ms after bounce. The large 
positive specific angular momentum gradient in its equatorial regions makes the left- 
hand side of eq. (HI) positive and convective overturn is confined to regions near the 
axis and large equatorial radii where the j-gradient flattens off. The shock in model 
s20.7r is deformed to oblate shape and its average radius is considerably greater than 
at the same time in the nonrotating model, but a top-bottom asymmetry typical for 
the SASI is absent. According to the results of Ott et al. [25|, rotation not only 
inhibits convection, but also delm^ and alters the growth of the axisymmetric SASI 
significantly. As pointed out by 3^, 31 1 , the situation may be different in 3D where a 
spiral-typc azimuthal SASI may develop in rapidly rotating SN cores. 

In nonrotating and slowly rotating core-collapse SNe, convection generically occurs in 
three distinct ways: (1) Prompt convection, driven by the negative entropy gradient 
left behind by the stalled shock and strongly dependent on the seed perturbations 
present in the core [1, H^l, (2) as lepton-gradient driven convection in the PNS (e.g.. 



12, Hi, [si), and (3) as the already mentioned neutrino-driven convection in the gain 



region behind the shock (e.g., 



2.2. GW Emission 

In nonrotating or slowly rotating core-collapse SNe and in the absence of strong 
PNS pulsations (discussed in section 2]), the dominant multi-D dynamics and, hence, 
the primary emission process of GWs is associated with convective overturn and the 
SASI. In addition, and also related to the asymmetry due to convection and SASI, 
anisotropic emission of neutrinos adds a slowly-varying component to the GW signal 
with characteristic frequencies of the order of 10 Hz. Recent quantitative estimates 
of the latter can be found in 0, [13, |3^. Here, I focus on the higher-frequency, more 
readily-detectable GW signal from fluid motions and present in the left-hand panel 
of figure [2] the GW signal from the postbounce simulations of Ott et al. [2^ (models 
s20.nr and s20.7r as discussed above) that were run with full 2D angle-dependent 
neutrino-radiation hydrodynamics. The GW signals include the contributions from 
PNS convection, neutrino-driven convection and SASI. Scaled to a fiducial galactic 
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Figure 2. Left: GW signal h+ (scaled by the source distance D and in units 
of cm) of convection in the PNS and neutrino-driven convection and SASI in the 
postshock region in the axisymmetric postbounce models s20.nr (nonrotating, 
black lines) and s20.7r (rapidly rotating, green lines) of [Ssj . Rapid rotation 
effectively damps convection and modifies the SASf, leading to a lower-amplitude 
and more slowly varying GW signal. Right: Characteristic strain spectra /ichar 
|36|| of the same GW signals assumed to be emitted at D = 10 kpc and contrasted 
with theoretical initial, enhanced, and advanced LIGO design noise curves. Note 
that, in practice, initial LIGO's low-frequency sensitivity is identical to that of 
enhanced LIGO [H. 



source distance of 10 kpc, the peak amplitudes correspond to a dimensionless GW 
strains of 1 .3 x 10~^^ and are reached by the nonrotating model s20.nr. Convection in 
the rapidly rotating model is strongly inhibited, both in the PNS and in the postshock 
region. This is reflected in its low-amplitude, slowly-varying GW signal. 

In the right panel of figure [21 I plot the characteristic gravitational wave strain 
of the two postbounce SN models, defined by Flanagan & Hughes [s^ as /ichar(/) = 
D^^ y/2TT^^Gc~^dEQ-\/vdf~^ , where D is the distance to the source (set to 10 kpc) and 
di?GW / df is the GW spectral energy density. In order to contrast the /ichar spectra 
with detector sensitivity, I also plot the design sensitivity h^ms of initial LIGO [s^, 
enhanced LIGO [3§|, and advanced LIGO (in burst mode) According to the 

simulation data of Ott et al., in a nonrotating or slowly- rotating SN, most of the 
energy in GWs is emitted in a broad frequency range of 300 — 1000 Hz (this is in rough 
agreement with by a combination of the high-/, small-scale PNS convection and 
downfiow plumes of material caused by the SASI that penetrate to small radii and 
are rapidly decelerated at the PNS surface [i , 24 1 . The emission due to convection in 
the rapidly rotating model is much weaker, lacks the high-/ component and occurs 
primarily at frequencies of 10 — 200 Hz. 

It is important to note two aspects: (1) The nature of the GW emission by 
convection and SASI is stochastic [4|,|35[ and the emitted GWs will be of uncorrelated, 
essentially random polarization. This is a consequence of the turbulent, chaotic 
emission dynamics that depend in a non-deterministic fashion on local conditions 
and can vary greatly from one SN to the next. (2) The detectability of the GWs from 
convection and SASI (i.e., the magnitude of /ichar in the right panel of figure[21) depends 
crucially on the duration of emission. The SASI and neutrino-driven convection shut 
off once an explosion has fully developed or a BH has formed. Hence, it may last 
~ 200 ms — 3 s, the upper limit depending on the onset of explosion, the accretion 
rate set by the progenitor and the EOS-dependent mass-limit for BH formation. PNS 
convection would be shut off by BH formation, but if an energetic explosion occurs. 
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it is likely to continue for many seconds until the PNS has reached neutrino-less /3- 
equilibrium. See 0] for scalings of the emitted GW energy with the duration of the 
convective signal. 



3. The Magnetorotational Mechanism 



Most (> 99%) core-collapse SN progenitors are expected to be rather slowly rotating 
and have central periods of tens to hundreds of seconds 4l|, |42|. However, there exists 
a variety of progenitor channels that could yield iron cores with precollapse periods 
of the order of seconds (e.g., 4^, 3). Conservation of angular momentum in collapse 
leads to a spin- up by a factor of ~ 1000 ji^]. Hence, an iron core with Pq ~ 2 s results 
in a PNS with P ^ 2 ms and a rotational energy of the rigidly-rotating PNS and of the 
strongly differentially rotating postshock layer of a few x lOB [42, 45|. A fraction of 
this energy would be sufficient to power an energetic core-collapse SN provided there 
is an efficient mediating mechanism for converting rotational energy into linear kinetic 
energy of the ejecta. 
Theory and early simulations [ii 



471 . |48| . |49| have shown that magnetorotational 



processes may constitute such a mechanism and can drive collimated outflows, leading 
to energetic jet-driven bipolar explosions [sot . Recently, this magnetorotational 
mechanism has received much attention (e.g., [U 



53l . |5j| and references 



51 



therein) owing to increasing observational evidence for collimated jets in the context 
of gamma-ray bursts (GRBs) and hypernovae (hyper-energetic > 10-B SNe) and the 
association of stellar collapse with the long-soft class of GRBs (e.g., jssf). 

Recent simulations [45, 52, HI, 54, 5^] suggest that magnetic fields of the order of 
10^^ G with strong toroidal components are required to yield the necessary magnetic 
stresses to drive a strong bipolar explosion. Rapidly rotating progenitors may have 
rather strongly magnetized cores, with toroidal fields of ~ 10^ — 10^^ G (the poloidal 
components are 1 — 2 orders of magnitude weaker) (4l] . l43l |. Flux compression in 
collapse scales cx (p/po)'^^^ S 111 and can lead to an amplification of both toroidal 
and poloidal components by up to a factor of 1000. Further amplification is provided 
for the toroidal com pon ent by rotational winding of poloidal into toroidal field (f2- 
dynamo, e.g., (isl. Isl-lsTt . a linear process ) an d by the exponentially-growing magneto- 
rotational instability (MRI, e.g., 51, 5^, 5^ and references therein) that has the 
potential to provide dynamically-relevant saturation field strengths of the order of 
10^^ G [H^. Both, fi-dynamo and MRI operate on rotational shear (dfl/dlnr) and 
convert energy stored in differential rotatior||| into magnetic field. Strong differential 
rotation in the outer SN core is a generic feature of rotating core collapse ji^, ■ 

The to-date most microphysically complete MHD core-collapse simulations (though 
in a Newtonian framework) were carried out by Burrows and collaborators 3, Sq . 
Under the provision that the MRI operates as expected 58[ and the postbounce 



rotational energy is of the order of ~ 10 B (implying Pq < 4 s), these authors 
find robust and powerful bipolar magnetorotational explosions that reach hypernova 
energies, leave behind a protomagnetar and could be the stage-setting precursor to a 
subsequent GRB 

II The lowest-energy state of a rotating fluid at fixed total angular momentum is uniform rotation. 
Hence, the difference in energy between uniform and differential rotation at fixed total angular 
momentum can be considered as "free energy of differential rotation." 
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Figure 3. Left: GW signal (rescaled by distance D and in units of cm) emitted 
by rotating core colla pse , bounce and postbounce nonaxisymmetric dynamics in 
model s20A2B4 of |60l . l6lll . Shown are the (top) and (bottom) polarizations 
as seen by observers situated in the equatorial plane (red lines) and along the 
polar axis (black lines). Note that the GWs are linearly polarized during the 
axisymmetric collapse and bounce phase and become predominantly elliptically 
polarized only tens of milliseconds after bounce. This panel is a variant of the 
right panel of figure 4 in [2. Right: Time-frequency (TF) analysis of the spectral 
energy density of the GWs emitted by model s20A2B4. The TF analysis was 
carried out with a 2 — ms Gaussian window that was shifted over the data in steps 
of 0.2 ms. Both the strong burst at core bounce and the energetic late postbounce 
emission leave clear and clearly separated marks in the TF diagram. 



3.1. GW Emission 

Without any reasonable doubt one can assume that rapid rotation is the primary multi- 
D component of the magnetorotational mechanism. Rapid rotation leads to a strong 
burst of GWs at core bounce, associated with the huge accelerations acting on the 
centrifugally deformed inner core. This GW emission process is the most extensively 
studied in the core-coUapsc context (see, e.g., Q for a review) and the currently 
best estimates of the GW signal come from 2D and 3D GR simulations that include 



a microphysical EOS and deleptonization during collapse [60|, |61|, |62, l63|. These 
studies have demonstrated that the GW signal from bounce is of generic morphology, 
exhibiting a small pre-bounce rise, a large negative peak roughly coincident with 
bounce and a subsequent ringdown signal emitted as the PNS settles to its new quasi- 
hydrostatic configuration. Comparing 2D and 3D simulations, 6^, 61 1 demonstrated 



that even extremely rapidly rotating cores stay axisymmetric through bounce, which 
constrains the GW emission to linear polarization (emission only in According 
to the results of [H^l , a core collapsing with a precollapse period Pq ~ 2 s {Pq 4 s) 
located at 10 kpc emits a peak GW amplitude of |/imax| ~ 1 x 10"^"'^ (|/imax| ^ 5x 10~^^) 
and a total energy in GWs of the order of a few x 10~^ Mqc^ with most of the emission 
being concentrated at 500 — 800 Hz, extending to lower frequencies with increasing 
rotation. 

Although realistic PNSs are unlikely to undergo the classical dynamical MacLaurin- 
type nonaxisymmetric instability at rotation rates > 27% [g^I, nonaxisym- 
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metric dynamics may still develop via a low-T/|Vl^| corotation instability [64|. For 
equilibrium NS models, this was first discovered by l65| and has since been found to 
occur also in more realistic postbounce SN settings [60l . 61, 67 1. 

In the left panel of figure [3l I plot the gravitational wave polarizations (top) and 
hx (bottom) as seen by equatorial (black lines) and polar observers (red lines) emitted 
by model s20A2B4 in the 3 + 1 GR framework of Ott et al. [63,[6lj. This model used 



a 20- Mq progenitor with am iron core set up to spin roughly uniformly at a period of 
~ 1 s. The purely axisymmetric (£ = 2, m = in terms of spherical harmonics) bounce 
signal is followed by a primarily axisymmetric ringdown. Nonaxisymmetric dynamics 
develops in the postbounce phase and becomes relevant some 40 ms after bounce as 
indicated by the rise of the GW signal emitted along the poles {£ — 2,m = 2) due 
to the quadrupole components of the nonaxisymmetric dynamics. The right panel of 
the same figure displays a time-frequency analysis of the GW spectral energy density. 
The quick change in the waveform at bounce leads to power in a broad range of 
frequencies, but most of the energy is emitted around 350—400 Hz in this rather rapidly 
rotating core. The nonaxsiymmetric component kicks in at higher frequencies around 
~ 900 — 950 Hz (with secondary, modulating components at ^ 700 Hz and ^ 500 Hz) 
and, despite emitting lower-amplitude GWs, is significantly more energetic than the 
bounce signal, emitting roughly i?GW ~ 1.6 x 10^ ''M0C^ until the simulation was 
terminated. The nonaxisymmetric dynamics could last longer, possibly for hundreds 
of milliseconds 0, 61 1 j but may be in competition with the MRI, an aspect that 
remains yet to be investigated. 



4. The Acoustic Mechanism 

The acoustic mechanism, proposed in a series of papers by Burrows and 
collaborators 0,0, (using 2D radiation- hydrodynamic simulations), is based on the 
excitation of PNS core pulsations (primarily £ = {1,2} 5-modes) by turbulence and 
by accretion downstreams through the unstable and highly-deformed stalled shock in 
the non-linear phase of the SASI. In the simulations of Burrows et al., fueled by the 
gravitational energy of anisotropic accretion, the PNS pulsations ramp up to non-linear 
amplitudes over many hundreds of milliseconds and damp via the emission of strong 
sound waves. These travel down the steep density gradient present in the postshock 
region, thus steepen to shocks and very efficiently deposit their energy behind and 
in the shock. Explosions in the acoustic mechanism are robust, but set in only after 
~ 1 s after bounce and tend to have energies on the lower side of what is observed. 

There are multiple caveats associated with the acoustic mechanism: (1) It is so- 
far unconfirmed, but also not yet ruled out by other SN groups (e.g., [2j|). (2) 
Perturbation theory suggests that the PNS mode amplitudes may be limited by a 
parametric instability involving high-order modes that are not presently resolved in 
numerical simulations [6^. (3) If the neutrino mechanism is effective, it will probably 
explode the star before the PNS pulsations have time to grow to large amplitudes. (4) 
As pointed out in section [21 rapid rotation has a stabilizing effect on convection and 
SASI and, hence, is likely to also inhibit the growth of PNS pulsations. In addition, 
rapidly rotating SNe are likely to explode by the magnetorotational mechanism much 
before PNS pulsations can reach large amplitudes. (5) It is not clear how the relaxation 
of axisymmetry would affect the power and spatial character of the PNS pulsations. 
3D simulations have yet to be performed. 
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Figure 4. Left: GW signals (rescaled by distance D and in units of cm) 
emitted in a set of representative nonrotating models from 0| that explode via the 
acoustic mechanism. At early times, the GW signal contains a burst associated 
with prompt convection at bounce and a longer-term broad-band low-amplitude 
contribution from PNS convection, neutrino-driven convection, and the SASI. At 
later times > 300 — 600 ms after bounce) , the GW emission becomes dominated 
by PNS core pulsations, leading to a narrowband signal (a zoomed-in view is 
shown in the inset plot). Around the time of the onset of explosion the largest 
GW amplitudes are reached. Right: Characteristic GW strain spectra /i^har 13611 
(at lOkpc) of the same models contrasted with initial, enhanced and advanced 
LIGO design noise curves. 



4-1- GW Emission 

The early postbounce GW emission associated with the acoustic mechanism is similar 
to that of the neutrino mechanism discussed in section 12.21 and comes primarily from 
convection, the SASI and SASI downflow plumes that are decelerated at the PNS 
core's outer edge. The latter are drivers of the PNS pulsations whose quadrupole 
components dominate the GW emission in the later postbounce phase. They are the 
most characteristic dynamical feature of the acoustic mechanism. 

In the left panel of figure [H I plot sample GW signals extracted from the simulations 
of Burrows et al. (a more detailed discussion of these models was presented in 
for a variety of progenitors. Present is a burst of prompt convection around core 
bounce, followed by a phase of neutrino-driven convection with little SASI. Around 
~ 300 — 600 ms after bounce, the non-linear SASI regime is reached, downflow plumes 
reach the PNS core and the core pulsations begin to grow and emit GWs at the 
frequencies of their quadrupole components (see inset plot) that arc set by the core's 
structure. While the GW emission in these models at intermediate times stems from 
£ — 2 harmonics of the then dominant £ = 1 mode, an £ = 2 eigenmode at higher 
frequencies appears dominant at times ^ 1 — 1.2 s, explaining the increase of the GW 
amplitudes around the onset of explosion. 

The PNS core pulsations of the acoustic mechanism are the strongest GW emission 
mechanism proposed in the stellar collapse context. Predicted strain amplitudes range 
up to h+ - 10-20 (at lOkpc, and total GW energies are 

£;gw - IO-^Mqc^ and 

above, emitted in a narrow frequency band set by the PNS pulsation frequencies 
and their temporal evolutions. The right panel of figure [51 contrasting /ichar spectra 
with detector design sensitivities, demonstrates that a galactic SN exploding via the 
acoustic mechanism would be difficult to miss even with first-generation detectors. 
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Figure 5. Scheme summarizing the associations of particular multi-D processes 
and characteristic GW emissions with the individual candidate core-collapse SN 
mechanisms. 

5. Discussion and Conclusions 

Each of the three SN mechanism that I consider in this article has an unique muIti-D 
dynamical process intimately linked with it. The acoustic mechanism, while involving 
neutrino-driven convection and SASI, is fundamentally reliant on non-linear PNS core 
pulsations to drive the acoustic SN engine. In the magnetorotational mechanism, on 
the other hand, rapid rotation is the key ingredient, damping convection and SASI, 
and facilitating B-field amplification and a bipolar jet-driven explosion. Finally, the 
neutrino mechanism's key multi-D aspect is convection and SASI. While PNS core 
pulsations could certainly be excited, they would not have time to grow to large 
amplitudes if the neutrino mechanism was effective and launched an early explosion. 

Now, as delineated in the previous sections and portrayed by figures HHH the three 
distinct multi-D processes dominant in the three SN mechanisms have very different 
GW emission characteristics: (1) Convection and SASI of the neutrino mechanism 
lead to a low-amplitude, broadband "stochastic" signal. (2) The rapid rotation 
necessary for the magnetorotational mechanism is reflected in GWs by a strong 
axisymmetric burst at core bounce and a possible prolonged narrowband postbounce 
GW signal from nonaxisymmetric dynamics. (3) The strong PNS core pulsations of 
the acoustic mechanism - if they obtain - emit copiously in narrow intervals about 
their characteristic frequencies. 

With the above information a simple picture can be drawn (figure [5]), relating 
a particular SN mechanism with a particular flavor of multi-D dynamics and an 
unique GW signature. Hence, with an observation of GWs from a core-collapse 
SN, the first strong observational constraints could be put on the SN mechanism. 
The magnetorotational and acoustic mechanisms should be clearly visible throughout 
the Milky Way with initial and enhanced interferometric detector technology while 
advanced and probably third-generation observatories (e.g., ET [t^I) niay be needed 
to constrain the SN mechanism for core-collapse events at 3 — 5 Mpc out to which the 
integrated SN rate is ^ 0.5yr~^ 0,[zH- Due to the low-amplitude, broadband nature 
of its GW emission, the neutrino mechanism is the most difficult to detect in GWs, 
but a non-detection of GWs from a galactic event with initial and enhanced detectors 
would quite clearly rule out the other two candidate mechanisms. 

The way in which I have portrayed the candidate SN mechanisms and their GW 
signatures in this article is idealized: I have focusscd only on the primary multi-D 
dynamics in each mechanism and its associated GW emission, neglecting contributions 
from secondary and, in particular, low-frequency components, such as GWs from 
anisotropic neutrino emission, magnetic stresses, or explosion asymmetries (see Q for 
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a detailed discussion). I have also neglected the possibility of the combined action of 
multiple mechanisms. For example, an explosion could be driven by a combination of 
neutrino heating and MHD processes in a moderately rapidly rotating SN. Yet, even 
in such a case, the dominant dynamical multi-D component would also dominate the 
GW emission and the simple scheme displayed in figure [5] is essentially preserved. 

Further theoretical and physically more complete and accurate computational work 
is necessary to gain a better understanding of the various possible SN mechanisms 
and of the details of their GW signatures. But, as 1 have attempted to show in this 
article, even with current knowledge, the GW signal of a galactic SN today may very 
likely provide strong hints for one mechanism and/or smoking-gun evidence against 
another. 
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